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ABSTRACT 

A rigorous analytical study of the dispersion relations of weakly amplified 
transverse fluctuations with wave vectors (k \\ B) parallel to the uniform back- 
ground magnetic field B in an anisotropic bi-Maxwellian magnetized electron- 
proton plasma is presented. A general analytical instability condition is derived 
that holds for different values of the electron (A e ) and proton (A p ) tempera- 
ture anisotropics. We determine the conditions for which the weakly amplified 
LH-handed polarized Alfven-proton-cyclotron and RH-handed polarized Alfven- 
Whistler-electron-cyclotron branches can be excited. For different regimes of the 
electron plasma frequeny phase speed w = u Pte /(kc) these branches reduce to 
the RH and LH polarized Alfven waves, RH polarized high- and low-phase speed 
Whistler, RH polarized proton and LH polarized electron cyclotron modes. An- 
alytic instability threshold conditions are derived in terms of the combined tem- 
perature anisotropy A = T±/T\\, the parallel plasma beta j3\\ = 8nn e kBT\\/B 2 and 
the electron plasma frequency phase speed w = w p , e / (kc) for each mode. 

The results of our instability study are applied to the observed solar wind 
magnetic turbulence at values of 90 < w < 330. According to the existence 
conditions of the different instabilities, only the left-handed and right-handed 
polarized Alfven wave instabilities can operate here. Besides the electron-proton 
mass ratio /i = 1836, the Alfvenic instability threshold conditions are controlled 
by the single observed plasma parameter w. The Alfvenic instability diagram 
explains well the main characteristic properties of the observed solar wind fluc- 
tuations. Especially, the observed confinement limits at small parallel plasma 
beta values are explained. 

Subject headings: plasmas - instabilities - turbulence - magnetic fields - solar wind 
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Introduction 



Kinetic plasma relaxation and turbulence generation processes are responsible for 
the observed properties of the solar wind plasma which is the only cosmic collisionpoor 



20091 1. Although the 



plasma accessible to detailed in-situ satellite observations (IBale et al. 
detailed plasma relaxation processes are not understood, the observed electron and proton 
distribution functions are close to bi-Maxwellian distributions with different temperatures 
along and perpendicular to the ordered magnetic field direction. Ten years of Wind/SWE 



data (IKasper et al. 



20021 ) have demonstrated that the proton and electro n temperature 



aniso tropics A = T±/T» are bounded by mirror and firehose instabilities ( IHellinger et al 
20061 ) at large values of the parallel plasma beta /3\\ = &imkBT\\/ B 2 > 1. In the parameter 



plane defined by the temperature anisotropy A = T±/T\\ and the parallel plasma beta (3\\, 
stable plasma configuration are only possible within a rhomb-like configuration around 
~ 1, whose limits are defined by the threshold conditions for the mirror and firehose 
instabilities. If a plasma would start with parameter values outside this rhomb-like 
configuration, it immediately would generate fluctuations via the mirror and firehose 
instabilities, which quickly relax the plasma distribution into the stable regime within the 



rhomb-configuration. The plasma pa rameters of other di 



interstellar and intracluste r medium (ISc 



compact, massive objects (ISharma et al. 



lekochihin et al. 



ute c osmic plasmas including the 



20051 ) and accretion disks around 



20071 ) are similar to the solar wind plasma, so that 



the temperature-anisitropy instabilities should also operate in these systems. 

The study of linear electromagnetic instabilities in a collisionless, homogoneous, 



magnetized, electron-proton plasma with temperature anisotropics has a 



review s we refer the intere st ed reader to Ch . 7 of the monograph by [Garyj ( 119931 ). 



(1992 


) and 


Marsch 


(2006 





ong 



hi story; for 



Hellinger et al. 



Cuperman 



( 120061 ) considered linear instability calculations 



for the oblique mirror, firehose and proton cyclotron instabilities, and represented the 
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apprq xi mate threshold conditi ons fr om the work of iGarv et al.l (119941 . l200ll ) , iGary fc Lee 



f ll994h 



Samsonov et al. 



feOOlD and 



Pokhotelov et al 



(2004) by analytic relations of the 



form A = 1 + a(/3|| — (3 ) b where a, b and (3 are fitted parameter values. 

In order to understand the confinement limits also at small values of the parallel 
plasma beta 0u < 1, we analyze here rigorously the full linear dispersion relation in a 
collisionless homogenous plasma with anisotropic (A ^ 1) bi-Maxwellian particle velocity 
distributions of electrons and protons for electromagnetic fluctuations with wave vectors 
(k x B = 0) parallel to the uniform background magnetic field B. With a typical solar 
wind temperature T = 10 5 T 5 K, the thermal particle energy k^T = 8.6T 5 eV is much less 
than the electron rest mass m e c = 5.11 • 10 5 eV, so that the use of the nonrelativistic 
linearized Vlasov/Maxwell equations is justified. For weakly (7 <C ur) amplified wave 
solutions we derive analytically existence and instability conditions, ujr and 7 here refer to 
real and imaginar y part of the com plex frquency u = + 17. Our analysis follows closely 



the recent study (jSchlickeiser 



2010|) for equal-mass pair plasmas - hereafter referred to as 



paper S. Our investigation is restricted to weakly amplified (7 <^ ujr) solutions covering 
the left-handed (LH) polarized Alfven-proton cyclotron branch and the right-handed 
(RH) polarized Alfven- Whistler-electron cyclotron branch. The corresponding analysis 
of weakly propagating (ojr 7) solutions, including mirror, firehose, electron cyctronic 
and cold magnetized Weibel fluctuations, is the subject of a subsequent paper. As we 
will demonstrate below, for the case of parallel (k x B = 0) wavevectors simple analytical 
threshold conditions for the two branches can be deerived in terms of the combined 
temperature anisotropy A = T±/T\\, the parallel plasma beta /3\\ = 8im e kBT\\/B 2 , the 
electron-proton mass ratio and the electron plasma frequency phase speed w = uj p ^/{kc). 
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2. Dispersion relations 



2.1. Basic equations 



For a nonzero background magnetic field strength the nonrelativistic dispersion 
relations for right-handed (RH) and left-handed (LH) pola rized fluctu ations with wave 



vectors k x B = in a thermal electron-proton plasma are (IGary 



19931 ) 



= D RHtLH (k,u) 



p.a 



a=p,e 



UJ 



_\/2ku at \\ 



= 0, 



(1) 



where we sum over a proton (p)-electron (e) plasma, and where k = \k\\\. 0J p ^ e denotes 
the electron plasma frequency, u a \\ = (kBT ay \\/m a ) l l 2 is the parallel thermal velocity of 
component a, Q a = e a B/{m a c) is the non-relativistic gyrofrequency, and A a = T aj ±/T aj \\ 
is the temperature anisotropy of component a, where the directional subscripts refer to 
directions relative to the background magnetic field. The dispersion relations ([T]) allow 
for different values of the proton and electron parallel temperatures and temperature 
anisotropies. 



Z(x) and Z'(x) denote the plasma dispersion function ( jFried fc Contdll96ll ) and its 
derivative 



Z(x) 



TT 



-1/2 



dt 



t — X 



(2) 



with the well-known properties 



Z\x) = -2[l + xZ(x)], 



(3) 



-6 - 



and 



Z{-x) = 2it 1 ' 2 ie- x2 - Z(x), Z'(-x) = Ait^ixe-* 2 + Z'{x) (4) 
We will frequently use the asymptotic expansions 

Z(x) ~ in 1/2 e~ x2 - 2x 

and 

Z(x) ~ W/V* 2 - - 

X 

where <r = if 3f(x) > 0, a = 1 if = and a = 2 if < 0. 

Paper S has demonstrated that the analytical analysis is enormously faoilitated if we 
work with space speed rather than frequencies. We therefore introduce the complex phase 
speeds 



2^ 
3 



\x\ -C 1 



(5) 



1 — — 

+ 2x 2 + 4a; 4 



b| » 1 



(6) 



/ = — = — : = R + ib, it=— , 6 = —, (7) 



the plasma frequency phase speed 



and the absolute value of the electron gyrofrequeny phase speed 

where \Q e \ — eB/m e c is the absolute value of the electron gyrofrequency. We also introduce 
the mass ratio 
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H = m p /m e = 1836, 
and the dimensionless proton and electron temperatures 



(10) 



2k B T e . 



m P c* 



1/2 



2k B T p 



1/2 



Swansonl (119891 ) we discuss high density 



Throughout this work, in the classification of 

plasmas with u Pie ^> \ fl e \, corresponding to w ^> b which applies to nearly all astrophysical 
plasmas. These plasmas are dense enough that the electron plasma frequency is much larger 
than the electron gyrofrequency, but small enough that elastic Coulomb collisions can be 
neglected. 

The two dispersion relations ([T]) then read 



w 
+ w z 



k 2 c 2 



f I /± 



a 



Lz ( fTb 



i /7± 

- (1 - Ap) Z' M 



-(i-A e) z' 



(12) 



We notice the symmetry A(—k\\,f) = A(k\\, f) of both dispersion relations allowing to 
consider only positive values of the wavenumber k > 0. In the following we will simplify 
the analysis by considering only equal parallel temperature plasmas (T e || = T^y) so that 



6 e = 9 and P = 0//i 



1/2 



The dispersion relations ffl2l can be separated into real and imaginary parts 
A(R, S) = &A(R, S) + i%A(R, S) = 0, implying the two conditions 



UA(R, S) = 0, %A(R, S)=0 



(13) 



-8 - 



In terms of the complex phase speed / = R + iS the real and imaginary parts of the two 
dispersion relations (fT2j) read 



= $tA RH>LH {R, S)=R 2 -S 2 -1+ 



w 



R m (fl 



1/2 







R + tS± 



+ 



2/i 



e 



R + tS± 



I'. 



4*z ( R+ t Tb ) + ( R+ t Tb 



w 2 S 



1 QZ^ 



1/2 



1/2 



e 



R + iS± - 



+ 



R + iSTb 




(14) 



and 



w 



= %A RH>LH (R,S) = 2RS+ — 



|V2 



, (R + iSTb 



R + iS±- 



w 2 S 







1 RZ'" 



1/2 



e 

1/2 



Here we consider solutions of the dispersion relations in the weak damping/amplification 
limit \S\ <C R- As shown in paper S in this limit the real part of the dispersion relation 
satisfies 



6 

1/2 



R + iS± 



R + iS ± - 



+ 3Z 



e 

R + iSTb 

e 



(15) 



9fcA(i2, 5 = 0) = 0, 
whereas the corresponding imaginary part is given by 



(16) 



S 



%A{R, S = 0) 

dftA(R,S=0) 
8R 



(17) 
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As before we introduce the parallel plasma beta 

2 w 2 8irn e kBT\\ 



0\\ 



b 2 B 2 ' 

which expresses the magnetic field strength as 



Gw 



1/2 



(19) 



in terms of the electron plasma frequency (w), the parallel temperature (G) and the parallel 
plasma beta /3m. Our restriction to high-density plasmas with u Pt£ > \Q e \ or w > b then 
requires to have parallel plasma betas (3\\ > 2 . 

3. Weakly damped and amplified solutions 

For weakly damped or amplified fluctuations Eqs. ( THj) - (fT5|) read 



o = s& RHtLH (R,s = o) = R 2 -i 



+ 



R nz ( ^ 



1/2 



e 



R± 



+ 



A (a 112 



2/i 







\ R ± b - 


)] 







and 



%A RH , LH (R,S = 0) 



w 



I- A 



1/2 



(20) 







R± 



+ 



w 2 R 







1/2 



1/2 







R± 



+ %Z 



-(l-A e )QZ' 
RT_b s 



,[R_±b 




(21) 



which have to be investigated for positive values of R > 0. 
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As noted before, these dispersion relations can be further reduced with the asymptotic 
expansions (J5j) and flBJ), depending on the absolute values of the arguments 



p±(R) 



1/2 



e 



S±(i2) 



|i?±6| 

e 



(22) 



of the plasma dispersion function and its derivative being small or large compared to unity. 
Scaling R = bx, corresponding to x = u)n/\£l e \, the absolute values of the arguments (j2"2"|) in 
terms of the parallel plasma beta read 



P±(x) 



Wfl 



1/2 



P 



1/2 



\x ± 



E±(x) 



w 



\x ± I 



(23) 



which are shown in Fig. 1 as a function of the normalized real frequency x. 



For parallel plasma beta values (3u {w 2 / ji) we note that P+(a;) ^> 1, E + (x) ^> 1 for 
all values of x, whereas P-(x) ^> 1 for x outside the small interval 



x 



A' 2 i 







1/2 



\L Wfl 1 / 2 ' fl Wfl 1 / 2 



(24) 



around the proton cyclotron frequency, and E-(x) 3> 1 for x outside the small interval 



x 



around the electron cyclotron frequency. 



w 



w 



(25) 



In the following we limit our analysis to such values of the parallel plasma beta 
<C (w 2 /[i), so that the asymptotic expansion ([6]) of the plasma dispersion function can 
be used except near the indicated proton and electron cyclotron frequencies. 

Irrespective of the values of P± and E± we note that both asymptotic expansions (jSJ) 
and fl6]) yield the same imaginary part of the dispersion relation 
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Fig. 1. — Absolute values of the arguments P±{x) and E±(x) of the plasma dispersion 
function ans its derivative in the weak damping/amplification limit as a function of x = R/b. 



= SArhMR, S = 0) = <n l/2 w 2 ^- 



+ 



A P 



± 1 



1 

=F - 



(R-b y2 

e © ; 



(26) 



4. Alfven, Whistler, cyclotron waves and electromagnetic light 

For parallel plasma beta values (3\\ <C (w 2 /fi) the asymptotic expansion (jBJ) yields 
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and 



R±b 

e 



e 



Ce 


VA 




v. 



e 



R±b 



e 2 



2( J R±6) 1 



(R±by 



e 



i + 



30 2 



2(R±b)' 



/iV2 (jR± A) 



1 + 



02 



b \2 



(27) 



(2f 



(29) 









0. 



02 



H{R± 



b\2 



1 + 



30 s 



2(i{R± 



b \2 



(30) 



We then obtain for the real part of the dispersion relations ( 1201 to lowest order in <C 1 
= &A rh ,lh{R, S = 0) = i? 2 - 1 - w 2 i? 



+ 



1 



RTb n(R±$ 



Q2 W 2 



+ 



1 - A 



_CR=F&) 2 (RTb) 3 \ 2/i 



2 w 2 



1 - Ar 



R 



With the scaling R = bx the dispersion relation f[3"Tj) reads 

(1 + /i)w 2 



= 3f?A WH (x,S = 0) 



6^ + 



2 



A e x ± (1 



-4, 



(1 ± /ix)(l =F a;) 
4p/ix =f (1 - Ap) 



x 2 - 1 



where we introduce the parallel plasma beta f fTSj) . 



^1 ± nx)' c 



(31) 



(32) 



In different limits the solutions of Eqs. f[3"Tj) and f l3"2"j) describe Alfven waves, Whistler 
waves, cyclotron waves and electromagnetic light. We consider each case in the next 
sections. 
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4.1. Alfven waves at phase speeds R <C for 2 <C f3\\ < w 2 / jj, 

For phase speeds R <^b/fi the dispersion relation (1311) simplifies to 



= ttk RH)LH (R,S = 0) = R 2 1 + 



2w 2 (l + /i) 

fo 2 



+4 

2 



1 

2R 



2-A e -A p ± =^(3 - A e ) =f ^(3 - A p ) 



i? 2 (l + ^)-[l + (A-l) 



(33) 



where we introduce the Alfven speed 



'1 + n)w 2 4irn e (m e + m p )c 2 c 



2 



b 2 B 2 V\ y 

the parallel plasma beta (I18p and the combined plasma temperature anisotropy 



(34) 



A = (35) 

For high-density plasmas Va <C c, the dispersion relation fl33|) yields the LH and RH 
polarized Alfven modes with the same 



which either can propagate forward and backward. Note that the condition R <C b/fi 
requires (w 2 /fi) ^> 1. 

Moreover, the four Alfven modes only exist for temperature anisotropics such that 
1 + (A — l)/3i| > corresponding to 
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A> 1 



1 

s 



(37) 



which includes the isotropic (A p = A e = A = 1) plasma temperature case. For small plasma 
betas < 1) the condition (j57|) is always fulfilled whereas for large plasma betas ((3\\ > 1) 
the combined anisotropy A has to be larger than 1 — (37 . 



Eq. (133|) also provides 



«9i? 



(38) 



so that for all four modes according to Eqs. (IT7|) and fl26l) the growth/damping rate is 



1/2&3 ( 1 



2(l + /i)9 J RV / i 1 /2 



Srh,lh 
R 



b [i 



i? 



± 1 



e e J 



(39) 



For isotropic (A p — A e — 1) plasma te mperatures al l four Alfven modes are damped in 
agreement with the general theorem of Irjrinca (199(3) on the electromagnetic stability of 



isotropic plasma populations. The Alfven damping rate in the isotropic case is given by 



Srh,lh(A p — A e — 1] 



e &z (n± v> +e — e- } 



2(1 + ^)9 L/i 1/2 

In order to drive the Alfven modes unstable, the growth rate (139!) has to be positive 
requiring that 



(40) 



r(A p ,A e ,P {] ) EE ± 

A 



'K/7. /.//-■ .,- • " ' — \ ^3/2 



- (A e - l)e 



+ X 



L// 



1/2 



e w n 



■(l±lix) 2 



+ A e e 



< 



This instability condition is analyzed further in the next section. 



(41) 
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4.2. LH-polarized Alfven-proton-cyclotron and RH-polarized 
Alfven- Whistler-electron cyclotron branches for small parallel plasma 

beta 0\i < 1 

Because we are primarily interested in low plasma beta plasmas we consider for 
completeness the solutions of the full dispersion relation (|32|) for subluminal solutions 
R <C 1 and small parallel plasma beta /3|| <C 1: 



(1 + fi)w 2 x 2 



o = m RH , LH (x,s = 0) 



1 ± (fj, — l)x — fix 2 



- 1 



yielding the quadratic equation 



(42) 



[(1 + fi)w 2 + fi\ x 2 =F (fx — l)x = 1 
Because \i = 1836 3> 1 this equation is well approximated by 



(43) 



i i 2\ 2 -,- 1 

1 + wja; =Fx = — 



with the solutions 



(44) 



XRH,LH 



covering the well-known (e.g. 



2(1 



Swanson 



in 2 ) 



1+ 4(l ± ^ ±1 



(45) 



19891 ) LH-polarized Alfven-proton-cyclotron and 



RH-polarized Alfven- Whistler-electron cyclotron branches for parallel propagation. 



For w <C ■>/ (/i/4) — 1 = 21.4, corresponding to large wavenumbers k 3> oo P)e /21Ac, we 
obtain 



1 1 

xrh ^ - + 



/i 1 + w 2 



(46) 
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and 



x LH ~ - 1 47 



The solution ( l47j) . subject to the constraint (124j) . represents the left-handed polarized 
proton-cyclotron waves. 

For values of w <C 21.4 the solution ( 146]) reduces to 



xrh ^ — — 2, (48 
1 + w z 

which for values of w <C 1 represents the right-handed polarized electron cyclotron waves 
subject to the constraint (1251) with the dispersion relation 



x RH -l-w 2 . (49) 
For intermediate values 1 <Cw <C 21.4 the solution f l48l) reduces to 



^ (50) 



W Z 



representing the right-handed polarized Whistler waves (R ~ b/w 2 ). 

For w ^> A/(/i/4) — 1 = 21.4, corresponding to small wavenumbers k <C w Pie /21.4c, we 
obtain for the solution fj45|) 

11 
^u(l + w 2 ) H 1/2 w 
which agrees with the Alfven wave solutions ( |36l) in the limit (3\\ <C 1. 

Eq. (T42|) for small plasma beta <C 1 also provides 
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d$tARH,LH(R, S = 0) (1 + /j,)w 2 x 



2±(fi-l)x 



2±(p-l)x 



OR 



(52) 



b [1 ± (/i — l)x — fix 2 } 2 6(1 + fi)w 2 x 3 ' 
We note that for all solutions ( 141)]) - f lBTl) . Eq. (j52J) is positive. According to Eqs. (IT7j) and 
(|26|) the growth/damping rate then is 



b{l + fi)w 2 x 3 
' 2±(n-l)x 



n 1 l 2 (l + fi)b 2 w A x 3 



Srh,lh 
1 



A 



X ± — 



1 

=F - 



9A M (-R,S = 0) 



e">' )2 + [i e (^l)±l]e- 



9[2±(/i-l)x] V/x 1 / 2 _ 
which generalizes the Alfvenic growth/damping rate ( 139]) to non-Alfvenic modes. Apart 



(53) 



Gary 



from t he different notation in phase speeds, the rate ( 153]) agrees with eq. (7.1.8) of 
(119931 ). For isotropic (A p = A e = 1) plasma temperatures all modes of the RH and LH 
branches, including the cyclotron and Whistler modes, are damped in agreement with 
Brinca's general theorem. We notice that, in order to drive the cyclotron and Whistler 
modes unstable, the growth rate (153]) has to be positive, requiring again the earlier derived 
instability condition fl4"T]) . 



4.3. Electromagnetic light at large frequencies R ^> b + 



For large frequencies R ^> b + the dispersion relation (j3T|) reduces to 



1 b w 2 b 2 

= ^A R h,lh(R^ = 0)^R 2 -1~(1 + -)w 2 tw 2 -- 



R R? 



~ R 2 - 1 - (1 + -)w 2 



1 



and 



, 1 & { A A p 



(54) 



o = qArhMR, s = o) * HJ I. < "" L " H 



V/^ 1/2 



(55) 
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so that for both polarisations to lowest order in (b/R) 2 we obtain the dispersion relation of 
electromagnetic light 



R 2 ~ 1 + (1 + -)w 2 



1 
/' 



with the same damping rate 



(56) 



_ tt^w 2 ( A p R 2 /@ 2 ~ R ye 2 

b ~ 2e~{j^ e +Aee 



(57) 



5. Instability conditions 

In order to drive RH-handed polarized Alfven-proton cylotron and the LH-handed 
polarized Alfven- Whistler-electron cyclotron branches unstable, the condition (1411) has to 
be fulfilled which, after multiplication with /i 3 / 4 , reads 



,'•1, 1 



with 



3/4 



e~ x p+Ay/ 4 e- X ' 



± ^ e -|ln M +ln(A p -l)-X2 



e |ln M +ln(A e -l)-X2 



(5? 



If U7 

= p-ji7l ( 1± ^)» X e = -^( 1 T^)- 

Both brackets in Eq. ( 1581) can be further reduced using the identities 



(59) 



e <ii+a2 _|_ g Q i— a 2 0^ a i 



2e ai cosh(a 2 ), e ai+a2 - e ai " a2 = 2e ai sinh(a 2 ), 



e c 1+C2 + e ci- C2 = 2e ci cosh ( C2 ) ; e 



ci+C2 gCi— a c o„ c i 



2e ci sinh(c 2 ) 



(60) 
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From the first bracket we identify 



a x = In [(A p - l)(A e - I)] 1 ' 2 - X2p + *\ , 



a 2 = —- In /i — In 



nl/2 



A. 



- W+ 



where 



W±(x) 



V2 Y~2 



W 



2 r 



±4x + U- l)(x 2 - -) 



From the second bracket we determine 



Cl = ln[A p A e f 2 -^±^., 



(-2 



- - In u — In 

4 



Ae 

A, 



1/2 



W 4 



With sinh(— 02) = — sinh(a2) and cosh(— C2) = cosh(c2) condition (158|) becomes 



±(A p -lf/ 2 (A e -l) 1 / 2 sinh 



- In a + In 

4 



A e -1 
A p -1 



1/2 



-A l J 2 A l J 2 x^ 2 C o^ 



- In u + In 

4 



.4 



1/2 



+ ^4 



(61) 



(62) 



(63) 



= ± (^ 3/4 (A e - 1) - fi- 3/ \A p - 1)) cosh W ± ± (fi 3/4 (A e - 1) + /^ 3/4 (A p - 1)) sinhiy ± 
- xfi 1/2 (// /4 A e + /i" 1/4 A p ) cosh W± - V /2 (/i 1/4 A e - /i" 1/4 A p ) sinh W± > 0, (64) 
which yields the general instability condition for the two branches in the form 



{A, - A p )[(^' 2 + 1) + (/i 3 / 2 - 1) tanh(W±] + (A + A P - 2)[(/i 3 / 2 - 1) + (/i 3 / 2 + 1) tanh(W ± ] 
: x/i ((A e + A p )[^/ 2 + 1 + (/iV2 _ l) tanh(W±)] + (A e - A p )[//V2 _ 1 + ^1/2 + 1) tanh(W±)]) 



> 1 



(65) 



Inserting the general dispersion relation for small plasma betas f|4"5l) for a; and W±(x) then 
provides the general instability conditions for the RH-handed polarized Alfven-proton 
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cylotron and the LH-handed polarized Alfven- Whistler-electron cyclotron branches, which 
to the best of our knowledge has not been derived before. 

For equal electron and proton temperature anisotropics (A p = A e = A ) the condition 
(|65|) reduces to 




- 1) + (/i 3/2 + 1) tanh(W±] 



> fix [/i 1/2 + 1 + (> 1/2 - 1) tanh(W±)] . (66) 
Note that for pair plasmas (// = 1) the condition ( 1661) reduces to Eq. (S-66) in paper S. 

Instead of working with the general dispersion relations (I45p used in the general 
instability condition (165 p . we will discuss different ranges of the plasma frequency phase 
speed w, where the dispersion relations (j4"5|) reduce to simpler limits, as demonstrated in 
Sect. 4.2. For w ^> 21.4 we can use the Alfvenic relation ( 1BT|) . whereas for small values 
of w < 21 A we can use the cyclotron and Whistler relations f|47|) . (115]). ( )49|) and (j50]l . 
respectively. 

Before investigating the individual modes, we inspect the instability condition (|66|) . The 
instability conditions holds for plasma betas (3\\ w 2 /fi, implying w 2 /2(3\\ ^> fi/2 = 918. 
Consequently, the argument of the tanh-function 



w 2 1 
7^f±(x), f ± (x) = ±4x + (f,-l)(x 2 -- 
2p\\ [i 



(67) 



is much larger than unity for f±(x) 3> 2j \x. 

The function f+{x) increases monotonically from its smallest negative value 
/+(0) = —(1 — = —0.9995 and becomes zero at 



x + 



/j+ 1 



1 - 



2/i 



1/2" 



0.95// 



-1/2 



-1/2 



(68) 
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Likewise, the function f-(x) attains its negative minimum value —1.0016 at xe 
2/(jU — 1) = 0.0011 and becomes zero at 



£ - 



1 + 



2/x 



1/2 



1.05// 



-1/2 



-1/2 



(69) 



(/i- 1)^/2 [ ^+1_ 

Hence only for values of x — R/b inside a small interval around /x" 1 ^ 2 , which for right-handed 
polarized fluctuations lies in the Whistler phase speed range, the functions f±(x) are smaller 
than 2/ ix. Both functions are well approximated by 



f±{x) ^ fo(x) = lix 2 - 1, 



(70) 



which is negative in the range of Alfven and proton cyclotron waves and positive in the 
range of electron cyclotron waves, respectively. 

5.1. Left-handed polarized Alfven waves and proton cyclotron waves 

In the range of Alfven and proton cyclotron waves x < n^ 1 the function f-(x) ~ — 1 is 
negative and practically constant, so that the instability condition (166 j) becomes 



A 



> XIL 



(ll 3/2 + 1) tanh 



w 



2 /3|| 



- (/. 3/2 - 1) 



(/i 1/2 + l)-(/i 1/2 -l)tanh 



w 



20|| 



(71) 



The tanh-function is well approximated by 



tanh(t) 

so that with 2/3||/w 2 < 2/n = 1/918 < 1 



1 + 



(72) 



-22 - 



tanh 



w 



i + 2 4 

1 ID* 



2A| 



we find for condition ( 17 ip 



1 



i - a + ^' 2 )4 



> fix 



l + (/i 1/2 -l)^ 



Because /i 3//2 3> 1 and /i 1 / 2 ^> 1 condition ( 1741) is well approximated by 



where 



V = — — — =/3\\//3c 



(73) 



(74) 



(75) 



(76) 



denotes a normalized parallel plasma beta value in terms of the critical plasma beta 



0c 



W 



(77) 



(1 + ^ 3 / 2 ) 

The restriction (3\\ <C w 2 / fi requires r\ <C fi 1 ^ 2 = 43. Obviously, we have to consider the two 
cases r] < 1 and 1 < rj < 43, corresponding to j3\\ < (3 C and (3 C < /3y < 43/3 c , respectively. 

For < j3 c the instability condition ([75]) becomes 



1 H + rj 
1 — > x- -. 

A 1-7]' 



(7i 



or 



1 + x)r/ 



(79) 



-23 - 



Likewise, for f3 c < /3m < 43/3 c we obtain 



corresponding to 



'— -l) >x t±l 

-a L)>x v -r 



(80) 



MP\\ > Pc) < 



V 



1: 



'1 + x)rj — (1 — /j,x) 

For plasma beta values equal to the critical value (rj = 1), the instability condition ( I75|) 
cannot be fulfilled. The two instability conditions ( !79|) and (IHTj) for the LH polarized Alfven 
and proton cyclotron waves are illustrated in Fig. 2 together with the instability conditions 
for the RH polarized Alfven and low phase speed Whistler waves, which we derive in the 
next subsection. 



We remark at this point that the threshold condition for the LH Alfven-proton 
cyclotron b ranch show n in F ig. 2 does not agree with the AIC threshold condition plotted 



in Fig. 1 of 



Bale et al. 



()2009|). 



5.2. Right-handed polarized Alfven waves and low phase speed Whistler waves 



In the range of the RH polarized Alfven and Whistler waves below x ji the 
function f+(x) ~ —1 is negative and practically constant, so that the instability condition 



together with 



-w 



! /2/3|| b 



ecomes 



An 



- 1 



(fi 3/2 + 1) tanh 



w 



20|l 



- (/. 3/2 - 1) 



> x\i 



(/i 1/2 + l)-(/i 1/2 -l)tanh 



w 



(82) 



With the approximation (1721) we derive for condition ( 182]) 
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_-l)(l- 77 )>x( M + 7 7 ) (83) 

^0 



For < (3 C the instability condition (1531 becomes 



^> < « < :?i-,)„ (84) 



Likewise, for j3 c < < 43/3 c we obtain 



For plasma beta values equal to the critical value (77 = 1), the instability condition ( |75|) 
cannot be fulfilled. The two instability conditions (184"]) and ( 185|) for the RH polarized Alfven 
and low phase speed Whistler waves are illustrated in Fig. 2. 



5.3. Right-handed polarized high phase speed Whistler waves and electron 

cyclotron waves 

In the range of high speed Whistler and electron cyclotron waves x ^T 1 ! 2 
the function f + (x) ~ fix 2 is positive. The instability condition (16"6"]) together with 
W + = w 2 f + {x)/{2(3\\) = fiw 2 x 2 /(2(3\\) then becomes 



1 



( M 3 / 2 -l) + (/x 3 / 2 + l)tanh 



2 2m 

pr 



2/?|| 



> 



(/x 1/2 + 1) + (/x 1/2 - 1) tanh 



2 2 
\iw x 



With the approximation ( 1721) we derive for condition ( 186]) 



(86) 



1 - 



1 - 



\IW 2 X 2 _ 



> X 



1 + 



Ai 



fXW 2 X 2 



*7) 
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Because x we find that 



so that 



Ai Ai -i 



1 . 1 + /xw4 a . n , 2 Al 



1 " X > * ^ ^ x[l + — (89) 

A i _ p n fiw z x z 



leading to 



* > rr^rr^i (90) 

We now discuss the instability conditions for individual wave modes. 



6. RH and LH polarized Alfven waves 

For w fi 1 / 2 > 21.4 we insert the Alfven wave dispersion relation ( f5Tj) in the 
conditions (!79|) . ( l8~Tj) . (|84|) and (|85|) . For < /3 C , corresponding to rj < 1, we obtain for the 
RH polarized Alfven waves 



A (a , R S^ w(l-7?) _ W(/3 C ~ Al) / Q1 \ 

and the LH polarized Alfven waves 



A ,o , a W -7/) _ ~ Al) / q9 N 

(W - 0^) - ( W + 0-1/2)^ ( w - ^1/2)^ _ ( w + 0-1/2)^,, ^ 

Likewise, for > /3 C , corresponding to r/ > 1, we obtain for the RH polarized Alfven waves 
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A (a ^ a W W(rj-1) _ w{P\\-Pc) , , 

and the LH polarized Alfven waves 



A (R u ^R)<- w(ry-l) _ 

°'^ lP|1 (W + H-V*)n ~(W- ^/2) ( w + ^-1/2)^,, _ ( w _ ^1/2)^ ^ 

These four Alfvenic instability conditions are shown in Fig. 3. Note, as indicated in 
the figure, that the limiting anisotropy values and parallel plasma betas for — > and 
—> oo and A — > oo, respectively, are solely determined by the mass ratio /i and the 
wavenumber-dependent plasma frequency phase speed w = co Pte /kc. 

The properties of the weakly amplified polarized LH and RH polarized Alfven mode 
are summarized in Table 1 and 2, respectively. 



7. Proton cyclotron waves 

For w <C 21.4 we insert the dispersion relation (j47D for proton cyclotron waves in the 
instability conditions f J79|) and (|8ip yielding 



A (R tt < R W ^-V) _ »(Pc-P\\) , . 



and 



MP]l Pc} < + ^ - (i + " (m + i)/3„ - (i + ^ 2 )/? c (%) 

These instability conditions cannot be fulfilled for anisotropics in the interval Aq G 
[/i/(l + /i), /i/(l + u> 2 )], including the isotropic case A = 1. The proton cyclotron instability 
conditions are illustrated in Fig. 4 as a function of the parallel plasma beta 
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Real phase speed range 


i? < 6/^ = u>9/(/i0j /2 ) 


Parallel plasma beta range 


e 2 < 0|| < w 2 //i 


Dispersion relation 


x = R/b= ^} 2w 


Existence condition 


w > = 43 > VW 4 ) - 1 = 21.4 


Instability conditions: 




w w 

w+fi- 1 / 2 ' UI— fl 1 / 2 


. _ [ 43u> ui 1 
Ll+43u>' U--43J 




A > 


w(/3c-/3||) 


for 0,, < C - (1+ - 3/2) 


(«,- M l/2) (3c _( tl , +At -l/2) ( g || 


for 0,, > C 


A ^ V>(P\\-tlc) 
^° ^ (« ) + M - 1 /2) ( g || _( u ,_ At i/2) ( g c 



Table 1: Properties of weakly amplified LH polarized Alfven wave mode 



Real phase speed range 


R < b/fi = we/(//0j /2 ) 


Parallel plasma beta range 


2 < 0|| < w 2 //i 


Dispersion relation 


x = R/b = -i- 


Existence condition 


w > /i 1 /^ = 43 > ^/(/i/4) - 1 = 21.4 


Instability conditions: 




U! W 
W+fl 1 / 2 ' U> — fl 1/2 


r ' u> 43u> "1 
Lu>+43> 43w-lJ 


for 0,1 < C = 7T f4m 




M/3c-/3||) 




(u.+ M 1 /2) (3c _( u ,_ M -i/2) ( g || 


for 0|| > C 


4 ^ M/3||-/3c) 
^° ^ (« J - M 1 /2) (3|r ( ly+M -i/2) ( g c 



Table 2: Properties of weakly amplified RH polarized Alfven wave mode 
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The properties of the weakly amplified proton cyclotron mode are summarized in Table 

3. 

8. Right-handed polarized low phase speed Whistler waves 

The Whistler-electron cyclotron dispersion relation ( 148]) yields values of x <C fi^ 1 ^ 2 
provided 



1 + w 



1 « (97) 



2 



which is equivalent to 



(/i 1/2 - 1) 1/2 = 6.5 < w < 21.4 (98) 
The instability conditions ( 1841) and ( |85l) then reduce to 



MPii < /3c) < = „ ,r: ^ 09) 



and 



(l + w 2 )(l 


- rf) 


(i+w 2 m- 




(1 + + w 2 ) 


— W 2 1] 


(l+fi + w 2 )/3 c 


- w 2 /3|| 


(1 + W 2 )(7] 


-1) 


(l + ^ 2 )(/3||- 





w z rj — (1 + fi + w 2 ) w 2 p\\ — (1 + /i + w 2 )p c 
These instability condition cannot be fulfilled for anisotropies in the interval 



1 + w 2 1 + w 2 , 

<A < ——, (101) 



1 + w 2 + jl w 2 
including the isotropic case Aq = 1. The RH low phase speed Whistler wave instability 
conditions are illustrated in Fig. 5 as a function of the parallel plasma beta f3\\. 
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The properties of the weakly amplified RH polarized low phase speed Whistler mode 
are summarized in Table 4. 



9. Right-handed polarized electron cyclotron waves and high phase speed 

Whistler waves 

The Whistler-electron cyclotron dispersion relation ( 1481) yields values of x ^ fx^ 1 ^ 2 
provided 



' > m~ 1/2 , (102) 



1 + w 2 
which is equivalent to 



w 



< (fi 1/2 - 1) 1/2 = 6.5 (103) 



Inserting the Whistler-electron cyclotron dispersion relation (14*8]) yields for the instability 
condition ( 190]) 



^w 2 {l + w 2 ) 
flW 4 — 2(1 + w z ) z p\\ 



which is illustrated in Fig. 6. 



The properties of the weakly amplified RH polarized high phase speed Whistler and 
electron cyclotron modes are summarized in Table 5. 
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Real phase speed range 




Parallel plasma beta range 


e 2 « A, « f 


Dispersion relation 


x-R/b- l [l 1+w2 } 


Existence condition 


w < VW 4 ) - 1 = 21A 


Instability conditions: 




1+fJ,' 1+w 2 






A > 


M(/?c-/3||) 


for 0|| < C - (1+ ; 3/2) 


(1+ W 2)/3 C -(M+1)/3|| 


for 0H > C 


(/i+ i )/3|| _ ( i + ^ )/3c 



Table 3: Properties of weakly amplified LH polarized proton cyclotron mode 



Real phase speed range 




Parallel plasma beta range 


2 < 0,1 < w 2 //i 


Dispersion relation 




Existence condition 


6.5 < w < 21.4 


Instability conditions: 






, ^ (l+«> 2 )(/3 c -/3||) 


for 0|| < C = (1+ - 3/2) 




for 0|| > C 


, ^ (l+«, 2 )(/3||-/3 c ) 
^0 > w 2 ^-(\+^ +w 2 )p c 



Table 4: Properties of weakly amplified RH polarized low phase speed Whistler mode 



Real phase speed range 


bfi" 1 / 2 < R < b 


Parallel plasma beta range 


9 2 < 0|| < w 2 /fi 


Dispersion relation 




Existence condition 


w < 6.5 


Instability condition: 




for 0|| < w 2 /fi 


A ^ nw 2 (l+w 2 ) 
^0 ^ M «,4_2(l+ W 2) 2/3|| 



Table 5: Properties of weakly amplified RH polarized high phase speed Whistler and electron 
cyclotron wave mode 
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10. Application to solar wind fluctuations 



The solar wind magnetic fluctuations measured by 
wavenumbers 



Bale et al 



(120091 ) near 1 AU have 



k ~ a/p p (105) 

with a = 0.56 ± 0.32 and the thermal proton gyroradius p p = 4.23 ■ 10 6 T 5 ly/2 B^ 1 cm, where 
we adopt an interplanetary magnetic field value B = 10~ 4 i?4 Gauss and a temperature 
T = 10 5 T5 K. With the solar wind particle density n e = 10 2 n2 cm -3 we find that the plasma 
frequency phase speed (E]) 



79.5 142 (T 5 n 2 )V 2 

W = ~ a B— = 1±0.57 B 4 (1 ° 6) 



covers the interval 



90^A^ — <w< 330^^ — (107) 

B4 -B4 



The electron gyrofrequency phase speed 



b = 3.13-10- 3 w-^ (108) 

is indeed much smaller than the electron plasma frequency phase speed w justifying the 
high density plasma approximation w ^> b. 

According to the existence conditions of the different instabilities for low plasma 
beta values summarized in Tables 1-5, only the left-handed and right-handed polarized 
Alfven wave instabilities can operate in the w- interval (I107p . Fig. 3 therefore applies to 
this application. Fig. 3 was calculated for a value of w = 280 lying within the range 
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(I107p of the observed solar wind turbulence, w = 280 yields for the critical plasma beta 
p c = w 2 /{l + fi^ 2 ) = 0.997. 

Besides the mass ratio \i the Alfvenic instability threshold diagram shown in Fig. 3 
is controlled by the single observed plasma parameter w. These two quantities define the 
characteristic plasma beta j3 c = w 2 / /i 3 / 2 and the parameter 

a = ~ =175- - — 109 

w — /I ' w 

As illustrated in Fig. 7 all asymptotic and characteristic values of the threshold anisotropies 
and the characteristic plasma beta values in the case of LH and RH polarized Alfven waves 
are solely determined by a and j3 c . 



11. Summary and conclusions 

We rigorously studied the dispersion relations of weakly amplified fluctuations with 
wave vectors k x Bo = in an anisotropic bi-Maxwellian magnetized proton-electron plasma 
by the appropriate Taylor expansion of the plasma dispersion function. Apparently for 
the first time, for equal parallel electron and proton temperatures, a general analytical 
instability condition ( 165]) is derived that holds for different values of the electron (A e ) 
and proton (A p ) temperature anisotropies. We determine the conditions for which the 
weakly amplified LH-handed polarized Alfven-proton-cyclotron and RH-handed polarized 
Alfven- Whistler-electron-cyclotron branches can be excited. For different regimes of the 
electron plasma frequeny phase speed w = uj P)B / (kc) these branches reduce to the RH and 
LH polarized Alfven waves, RH polarized high- and low-phase speed Whistler, RH polarized 
proton and LH polarized electron cyclotron modes. The properties of the individual 
modes are summarized in Tables 1-5. Analytic instability threshold conditions are derived 
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in terms of the combined temperature anisotropy A = T±/T», the parallel plasma beta 
(3\\ = 87m e k B T\\/ B 2 and the electron plasma frequency phase speed w = u p>e /(kc) for each 
mode. 

For large values of w ^> 21.4, corresponding to small wavenumbers k <C u p , e /2lAc, RH 
and LH polarized Alfven waves are excited in different regions of the temperature anisotropy 
(A)-versus parallel plasma beta (/3||)-parameter plane. For appropriate temperature 
anisotropies, RH polarized low and high phase speed Whistler waves are excited for 
intermediate values of 6.5 Cw< 21.4 and lCt«< 6.5, respectively. At values w <C 21.4 
LH polarized proton cyclotron waves can be excited, whereas RH polarized electron 
cyclotron waves can be e xcited at small values of w < 1. In agreement with the general 



theorem of 



Brinca 



(1l990l ) on the electromagnetic stability of isotropic plasma populations 



none of these modes can be excited for isotropic plasma distributions (A e — A p — 1). 



We ap ply the results o 



turbulence (IBale et al. 



our instability study to the observed solar wind magnetic 
20091 ). corresponding to values of 90 < w < 330. According to the 
existence conditions of the different instabilities, only the left-handed and right-handed 
polarized Alfven wave instabilities can operate here. Besides the electron-proton mass ratio 
fi = 1836, the Alfvenic instability threshold conditions are controlled by the single observed 
plasma parameter w. All asymptotic and characteristic values of the threshold anisotropies 
and the characteristic plasma beta values are solely determined by the values of a = 43 /w 



and (5 C 



w 



''//j?/ 2 . Comparing the Alfveni c 



observations shown in Fig. 1 of 



Bale et al 



insta bility diagram shown in Fig. 7 with the 



(120091 ). one notices that the main characteristic 



properties of the observed solar wind fluctuations are well reproduced by the instability 
conditions of parallel propagating Alfven waves. Especially, the observed confinement limits 
at small parallel plasma beta values are explained. It remains to be investigated in future 
work how weakly propagating and/or obliquely propagating weakly amplified modes add to 
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this diagram. 

We thank the referee for his constructive comments that helped to improve the 
manuscript. This work was supported by the Deutsche Forschungsgemeinschaft through 
grants Schl 201/19-1 and Schl 201/21-1. 
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Fig. 2. — Anisotropy diagram for LH polarized Alfven and proton cyclotron waves and for 
RH polarized Alfven and low phase speed Whistler waves for the case of equal electron and 
proton temperature anisotropies A e = A p = A . 77 = / (w 2 / fi) denotes the normalized 
plasma beta. Unstable regions are marked by "u", 
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Fig. 3. — Anisotropy diagram for LH and RH polarized Alfven waves for the case of equal 
electron and proton temperature anisotropies A e = A p = A . Stable regions are marked by 
" S" , unstable regions are marked by " u" . A value of w = 280 is adopted so that the results 
hold for parallel plasma beta values (3\\ < 5.4. 



10 3 



10 2 



H/Cl+w 2 ) 
10 1 




M/(l+n) 



icrp c 10- 



Fig. 4. — Anisotropy diagram for LH polarized proton cyclotron waves for the case of equal 
electron and proton temperature anisotropies A e = A p = A . Unstable regions are marked 
by "u", A value of w = 10 is adopted so that the results hold for parallel plasma beta values 
P\\ < 0.054. 
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Fig. 5. — Anisotropy diagram for RH polarized low phase speed Whistler waves for the case 
of equal electron and proton temperature anisotropies A e = A p = A . Unstable regions are 
marked by "u". A value of w = 10 is adopted so that the results hold for parallel plasma 
beta values (3\\ < 0.054. 
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Ao 



1+w' 




Fig. 6. — Anisotropy diagram for RH polarized high phase speed Whistler and electron 
cyclotron waves for the case of equal electron and proton temperature anisotropies A e = 
A p = A . Unstable regions are marked by "u", Values of w = 3 and w = 0.01 are adopted 
for high phase speed Whistlers and electron cyclotron waves, respectively, so that the results 
hold for parallel plasma beta values below 5 • 10 -3 and 5 • 10~ 6 , respectively. 
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Fig. 7. — Anisotropy diagram for LH and RH polarized Alfven waves in terms of the plasma 
parameters f3 c and a = jj 1 ^ 2 /w. Unstable regions are matrked by "u", 
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